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INTRODUCTION 


This document is the second report concerning the Ohio Land 
Allocation Model. During Phase I, reviews were made of models 
of land use change and the uses of LANDSAT data for modeling.^ 

A standard set of land categories was formulated and an initial 
set of tax assessment and tax parcel models were derived. 

There were several additional tasks carried out in this 
phase of the study. These include: 

1) The formulation of tax assessment models for the 
agricultural sector; 

2} Sampling of parcel size data from three central 
Ohio counties to establish the feasibility of 
converting the tax parcel data to land acreage 
data; 

3) Interpretation of aerial photogarphy for Franklin 
County, Ohio, for use in quantifying the potential 
errors in LANDSAT data interpretations; 

4) The formulation of a standardized data base for 
use in Phase III - checking the errors in LANDSAT 
interpretations , 

The results of these tasks are sxunmarized in this report. Also 
included are two appendices. Appendix A discusses the evaluation 
of the DEMOS model while Appendix B sxammarizes a study of a step 
down population projection method for use in Ohio. 


^Oscar Pisch and Steven I. Gordon, The Ohio Land Allocation Model ; 
Report on Phase I f Columbus, Ohio: Ohio Department of Economic 

and Community Development, January, 1976) . 
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Finally', the research design for the third and final phase 
of the land allocation study is given. This section describes 
how LANDSAT data interpretation accuracy will be tested against 
the aerial photograph interpretations for Franklin County. 


Chapter I 


Land Use and Tax Models for the Agricultural Sector 

Phase I of the Ohio Land Allocation Model established a number 
of cross-sectional statistical models of tax assessment and tax 
parcel changes in the residential, commercial, and industrial land 
use sectors. Each of these equations utilized the assessed value 
or number of parcels in each land category as the dependent vari- 
ables and employment and population as the independent variables. 
Equations were derived which explained tax base in a static sense 
Ci.e. 1967 residential assessed value for 88 Ohio counties as a 
function of 1967 employment and population) and in a dynamic sense 
(change in tax assessed value 1967 to 1972 as a function of change 
in employment and population during the same period) . 

Similar models for the agricultural sector did not work as 
well as those for the other land use sectors. This is because 
population and employment variables do not incorporate the factors 
important for agricultural production. For this reason, a data 
base relating to agricultural production was compiled. Data on 
cash receipts to Ohio farmers were available from the Ohio Agricul- 
tural Research and Development Center.^ Data on the acreage 
harvested for each crop type were extracted from Ohio Agricultural 
Statistics while data on the number of farms, average acreage 

^ Ohio Farm Income, 1967, 1972 . Ohio Agricultural Research and 
Development Center, Wooster, Ohio. 
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and land in farms by county for 1964-1974 were made available 
by the Ohio Crop Reporting Service.^ Table 1 lists the variables 
which were used in the analysis of the agricultural sector. 

Data Analysis 

The first step in the project was to analyze the data base 
using descriptive statistics. The results of this analysis give 
an overview of the general trends in Ohio's agricultural sector. 

Land use in Ohio is diversified. Of the state's 26.2 million acres 
of land, 11 percent was dedicated to urban and built up areas in 
1970; '65.3 percent was in farmland.; and 23.7 percent was cate- 
gorized as other land not in farms. ^ As evidenced from this data, 
the agricultural sector is clearly a dominant land use within 
the State of Ohio. 

Changes in the state's agriculture indicate uhat Ohio is fol- 
lowing the national trend in that the number of farms and land in 
farms is decreasing, while average acreage is increasing. Although 
Ohio is following the national trend, data indicate that the state 
is not doing so in as radical and rapid a pace as the national 
average. For example, the average national farm size in 1973 was about 
385 acres while Ohio's average farm size for the same year was 

^ — 

Ohio Crop Reporting Service, Ohio Agricultural Statistics, 1967 , 

1972 , Ohio Crop Reporting Service, unpublished data. 

4 

Ohio Farm Bureau Federation, "Land Use Task Report," December 4, 

1974. 

s * * * 

U.S. Department of Agriculture, Statistical Reporting and Economic 
Research Service, 1973. 


Table 1 


LIST OF VARIABLES. AGRICULTURAL SECTOR TAX BASE EQUATIONS 


V67001 - Cash Reccilpts 000' s $ Total Livestock 1967 

V67002 - Cash Receipts 000 's $ Dairy 1967 

V67003 - Cash Receipts 000' s $ Cattle 1967 

V67004 - Cash Receipts 000 's $ Hogs 1967 

V67005 - Cash Receipts 000' s $ Poultry 1967 

V67006 - Cash Receipts 000' s $ Sheep 1967 

V67007 - Cash Receipts 000 's $ Other 1967 

V67008 - Cash Receipts 000 's $ Total Crops 1967 

V670C9 - Cash Receipts 000' s $ Corn 1967 

V67010 - Cash Receipts 000 's $ Soybean 1967 

V67011 - Cash Receipts 000' s $ Wheat 1967 

V67012 - Cash Receipts 000' a $ Oats and Hay 1967 

V67013 - Cash Receipts 000' s $ Greenhouse 1967 

V67014 - Cash Receipts 000' s $ Vegs. 1967 

V67015 - Cash Receipts 000 's $ Other Crops 1967 

V67016 - Acres Harvest Coim For Grain 1967 

V67017 - Acres Harvest Soybeans for Beans 1967 

V67018 - Acres Harvest All Wheat 1967 

V67019 - Acres Harvest Oats for Grain 1967 

V67020 - Acres Harvest All Hay 1967 

V67021 - All Cattle and Calves (Head) 1967 

V67022 - All Hogs and Pigs (Head) 1967 

V67023 - All Sheep (Head) 1967 

V67024 - All Milk Cows (Head) 1967 

S67 Acres - Total Acres Harvested 1967 

LF 67 - Land in Farms 000 *s Acres 1967 

NF 67 - Number of Farms 1967 
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Table 1 (cont'd) 


AA 67 - 

INDO 20 

IIIDO 60 

V72D01 

V72002 

V72003 

V72004 

V72005 

V72006 

V72007 

V72008 

V72009 

V72010 

V72011 

V72012 

V72013 

V72014 

V72015 

V72016 

V72017 

V72018 

V72019 

V72020 

V72021 

V72022 

V72023 

V72024 


LIgy OF VARIABLES. AGRICULTURAL SECTOR TAX BASE EQUATIONS 


Average Acreage 1967 
- Aseeased Value agricultural Land 1967 
~ 'Xotal Assessed Value Agricultural Land 1967 
• Cach l^.ecelpts 000' a $ Total Livestock 1972 
Cash ?:i.cel;.)te 000* s $ Dairy 1972 
Cash Recelrta 000 *s $ Cattle 1972 
Cash Receipts 000 *s $ Hogs 1972 
Cash Receipts 000' s $Poultry 1972 
Cash Receipts 000' s $ Sheep 1972 
Cash Receipts 000' a $ Other 1972 
Cash Receipts 000 's $ Total Crops 1972 
Cash Receipts 000* s $ Corn 1972 
Cash Receipts 000' s $ Soybean 1972 
Cash Receipts 000 *s $ Wheat 1972 
Cash Receipts 000 's $0ats and Hay 1972 
Cash Receipts 000' s $ Greenhouse 1972 
Cash Receipts 000' s $ Vegs. 1972 
Cash Receipts 000 's $ Other Crops 1972 
Acres Harvest Corn For Grain 1972 
Acres Harvest Soybeans For Beans 1972 
Acres Harvest All Wheat 1972 
Acres Harvest Oats For Grain 1972 
Acres Harvest All Hay 1972 
All Cattle and Calves (Head) 1972 
All Hogs and Pigs (Head) 1972 
All Sheep 1972 (Head) 

All Milk Cows 1972 (Head) 
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Table 1 (cont'd) 

LIST OF VARIABLES, AGRICULTURAL SECTOR TAX BASE EQUATIONS 


S72 Acres 
LP 72 
NF 72 
AA 72 
INDO 25 
INDO 65 
SMSA 67 
SMSA 72 
AVLIV67 
AVCRP67 
AVLIV72 
AVCRP 72 
1021 
1022 
1023 
1020 
1060 
CHGCRP 
CHGLIV 


- Total Acres Harvested 1972 

- Land In Farms 000 *s Acres 1972 
~ Number of Farms 1972 

- Average Acreage 1972 

- Assessed Value Agricultural Land 1972 

- Total Assessed Value Agricultural Land 1972 

- Counties in an SMSA, 1967 

- Counties in an SMSA, 1972 

- Mean of Cash Receipts Livestock, 1962-1967, 000 's $ 

- Mean of Cash Receipts Crops, 1962-1967, 000 *s $ 

- Mean of Cash Receipts Livestock, 1968-1972, 000 *s $ 

- Mean of Cash Receipts Crops, 1968-1972, 000 's $ 

- Delta ind. Land Assessed Value, 1967-1972 

- Delta comm. Land Assessed Value, 1967-1972 

- Delta res. Land Assessed Value, 1967-1972 

- Delta ag. Land Assessed Value, 1967-1972 

- Delta & Total ag. Assessed Value, 1967-1972 

- AVCRP 7 2 minus AVCRP 6 7 

- AVLIV72 minus AVLIV67 
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almost 149 acres. From 1964-1974 the average acreage on Ohio farms 
increased by 7.2 acres. The 1969 Census of Agriculture reports 
that 96.2 percent of all farms in Ohio are less than 499 acres. 

Small farms, 1 to 99 acres constitute 46.5 percent of this total.® 
This data would seem to indicate that although Ohio is slowly fol- 
lowing the national trend, the state remains dominated by relatively 
small farms and small farmers as opposed to large corporate farming 
operations that are prevalent in other states. Tables 2 and 3 
illustrate these data. 


Table 2. Ohio 1964, 1969, 1974 


Land in Farms 
OOO's 

1964 1969 1974 

18,145 17,700 17,400 


Average Acres per 
Farm 

1964 1969 1974 

141.6 148.1 148.8 


Ntimber of Farms 
000 's 

1964 1969 1974 

131 120 117 


Source: Ohio Crop Reporting Service 


Table 3. United States 

Average Acres per 
Farm 

1964 1969 1974 

332 369 385 

Source: U.S. Department of Agriculture, Statistical Reporting 

Service, 1973. 


U.S. Department of Commerce, Bureau of the Census, Census of 
Agriculture, Ohio 1969. 


- 8 - 


7 ? 



The assessed value of farms in Ohio has increased 23.1 percent from 
1964 to 1969. Crops increased 15.4 percent and livestock, 

7 

poultry, and related products increased 29 percent. Data for 
1967, one of the years that will be used in this study, show that 
livestock products are 57 . 5 percent of the total while crops con- 

Q 

stitute 42.5 percent of the total agricultural production in Ohio. 

In 1972, the other year relevant to this study, livestock pro- 
ducts constituted 52.8 percent and crops proved to be 47.2 per- 
cent of the total. ^ Table 4 illustrates livestock and major 
crops by type. 

Table 4 



1967 

1972 

Total Livestock 

57.5% 

52.8% 

Cattle & Calves 

14.8% 

16.6% 

Hogs 

14.5% 

13.4% 

Dairy Products 

18.8% 

16.9% 

Total Crops 

42.5% 

47.2% 

Corn 

11.4% 

11.9% 

Soybeans 

12.3% 

18.5% 

Wheat 

5.3% 

3,7% 

1 

Census of Agriculture , 

1969. 



®Ohio Agricultural Research and Development Center, Ohio Farm 
Income 1967. 

9 . 

Ohio Agricultural Research and Development Center, Ohio Parra 

Income 1972. 
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The next step in the analysis was to delineate the inter- 
relationships among the variables. This was critical in that 
strong correlations among independent variables in the regres- 
sion equations (multicollinearity) can yield erroneous results. 

Several sets of variables were found to be correlated. 

First, the analysis showed the relationship among the cash 
receipts variables and the acres harvested variables to be col- 
linear. This mecins that in almost all instances the cash receipts 
from a specific crop and the acres harvested from that same crop 
will correlate significantly. Table 5 illustrates the correla- 
tion coefficients for cash receipts and acres harvested; For 
example, cash receipts for wheat (V72011) and acres harvested 
for wheat (V72018) show a correlation coefficient of .90505.' 
Similarly, cash receipts for soybeans (V72010) and acres harvested 
for soybeans (V72017) correlate at .99385. As a result of this 
collinearity the regression statements which included both cash 
receipts and acreage for the same crop would give biased results. 
It is for this reason that the regressions use only the cash 
receipts for each crop and livestock variable. The decision to 
eliminate the acres harvested variables was based on the fact 
that cash receipts proved to be a stronger explanation of assessed 
value . 
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Table 5 


V67009 

V67010 

V67011 

V67012 

V67016 

V72009 

V72010 

V72011 

V72016 

V72021 


V67018 
. 90000 
.85954 
.98199 
.65113 


RELATIONSHIPS AMONG CASH RECEIPTS AND ACREAGE 
HARVESTED VARIABLES 
CORRELATION COEFFICIENTS* 

V67019 V72QQ3 V72Q16 V72Q17 

.55295 

.55772 

. 63133 

.80253 


V72018 


.89848 .87421 .80287 

.84860 .99385 .84102 

.86372 .93204 .90505 

.67246 

.82175 


♦Underlined coefficients indicate variables which are collinear 
and should not be used in the same regression equation. 


Collinearity was also a factor in determining the relation- 
ship among major crop types. This is shown in Table 6. The 
analysis showed that certain types of crops were intimately 
related. One explanation for this phenomenon is that certain 
crops are grown together or alternatively in a rotation. Some 
crops require the same type of storage soil^ or involve similar 
transportation costs. Soybeans (V72010) and wheat (V72011) cor- 
relate significantly. Corn (V72009) is collinear with soybeans 
(V72010) and with wheat (V72011) . Based on the theory of trans- 
portation costs one would not expect vegetables (V72014) to 
correlate highly with wheat or corn. This is, in fact, the 
case. It may also be posited that vegetables are intensive while 
wheat and corn are extensive crops, hence they would not be 
highly related in terms of where they are grown. 

Given adjustments in the input data to avoid the use of 
collinear variables, analysis was carried out using linear regres- 
sion. The results are summarized in the next section. 

Results of Agricultural Sector Models 

The goal of the regression analysis was to create predictive 
models of the tax base related to production in the agricultural 
sector and the conversion of land to urban uses. In each equation, 
one of the agricultural tax base variables (e.g, IND060 - total 


of THb 
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Table 6 


V67010 

V670U 

V72010 

V72011 


RELATIONSHIPS AMONG MAJOR CROP TYPES 
CORRELATION COEFFICIENTS 


V67009 

.8918S 

.91392 


V67011 

.89100 

1.00000 


V72009 

0.89816 

0.88707 

.90784 

.90587 


V72011 

0.93534 

0.95602 

.94601 

1.00000 
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agricultural assessed value, 1967) , is the dependent variable 
and production variables or assessment variables in other 
categories (residential, commercial, industrial) are the in- 
dependent or explanatorv* variables. 

Tables 7 to 12 illustrate the results of these analyses. 

All the tables are similarly organized. Using Table 7 
as an example, one can see how to interpret the results. This 
table represents the models for assessed value of agricultural 
land as given by the title. Results of two equations are shown 
in this table - one for 1967 as shown in the left half, and one 
for 1972 as shown in the right half. On the left margin are 
the names of all the independent variables for the 1967 equa- 
tion. The first numbers following this are the B coefficients, 
the numbers by which one multiplies to arrive at the predicted 
value &r the dependent variable. In parentheses after this 
number are the t ratios which show the statistical significance 
of each of the B values. At the bottom of the table are shown 
the constant for each equation and the coefficient of determina- 
tion (R^) , The R^ and R^ (a) (adjusted) values show the propor- 
tion of the variance explained by the equation. 

Going back to Table 7 , we can illustrate how one can trans- 
late into a predictive equation. Using the 1967 equation: 
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IND025 » .921 (V67004) + 3.789 (V67007) + 1.718 (V67009) 

+ 1.008 (V67010) + 10.493 CV67012) + 1.177 CV67013) 
+ 2927.345 (SMSA67) + 3361.156 

If one wishes to predict the assessed value of agricultural 
land for any county in Ohio using this equation, one need only 
have the production figures for each of the crop types in the 
equation and know whether or not the county was in an SMSA. One 
can interpret each of the tables in a similar fashion. How 
these equations will be translated into models for use on the 
state computer system is discussed further in the final chapter 
of this report. 

Several interesting trends are shown by these regression 
equations. Tables 7 and 8 show that a strong relationship 

exists between agricultural production and assessed value in 

2 

both 1967 and 1972. The R (a) values for these equations 
range from .683 to .809. Certain products are shown to be more 
important in Ohio agriculture. These include hogs, corn, 
soybeans, greenhouse, and other crops categories. ^ 

Viewing the coefficients in these equations, it appears 
that there were changes in the importance of particular crops 
between the two years. Thus, the b coefficient for soybeans 
in 1967 is 1.008 while it is 1.808 in 1972. It may not be true, 
however, that these two single years are representative of long 


Table 7 


Assessed Value Agricultural Land* 

1967 1972 


Indep .Variables 

B 

t ratio 

B 

t ratio 

Indep. Var 

V67004 

.921 

(3.51) 

1.524 

(3.95) 

V72003 

V67007 

3.789 

(1.44) 

29.400 

(4.55) 

V72007 

V67009 

1.718 

(3.14) 

1.808 

(8.96) 

V72010 

V67010 

1.008 

(1.86) 

1.539 

(3.03) 

V72013 

V67012 

10.493 

(3.97) 

5003.686 

(3.41) 

SMSA72 

V67013 

1.177 

(3.35) 




SMSA67 

2927.345 

(3.03) 




C (constant) 

3361.156 


2952.364 


C 


.826 


.705 




.809 


.691 


R^ (a) 

Mean 14545.002 


19894. 

229 


A Mean 


5349.227 





* All coefficients are significant at the. 001 level except 

V67007 (.2) 

V67010 (.1) 
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Table 8 


Total Agricultural Assessed Value* 
1967 


1972 


Indep . 

Variables B 

t ratio 

B 

t ratio 

Indep. Var 

V67007 

19.228 

(4.21) 

1.102 

(2.79) 

V72002 

V67009 

3.508 

(7.62) 

1.873 

(3., 11) 

V72003 

V67012 

18.474 

(4.14) 

32.738 

(3.30) 

V72007 

V67013 

1.572 

(2.50) 

1.890 

(5.36) 

V72010 




7.146 

(1.30) 

V72012 




2.205 

(3.15) 

V72013 

SMSA67 

3827.692 

(2.21) 

6387.938 

(3.07) 

SMSA72 

C 

8781.813 


5405.612 


C 

r2 

.729 


.705 



R^(a) 

.716 


.683 


R^(a) 

Mean 

24880.412 


30561.794 


Mean 

A Mean 


5681.38 





*All coefficients are significant at the .01 level. 

Except V67013 (.02) 

SMSA67 (.05) 

V72012 (.2) 
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term trends in Ohio agriculture. For this reason data on the 
axmual cash receipts for crops and livestock were collected 
for 1962 through 1972. Thus^ instead of employing the single 
year cash receipts variable, a six and five year average (1962-* 
1967, 1968-1972) were input into the regression equations. 

The results are shown in Tables 9 and 10. Here, one can com- 
pare the results using the simple totals for one year (V67001, 
V67008) to the results using the averages (AVLIV67, AVCRP67). 
Not only do the b coefficients change but also they become 
stable over time. The coefficient for AVLIV67 is .632 and 
for AVLIV72 .648 - not significantly different. What this 
means is that changes in agricultural assessed value from 1967- 
1972 are accounted for almost entirely by the dummy variables 
for urbanization - SMSA67 , SMSA72, In other words, competition 
from urban land uses in SMSA counties was bidding up the price 
and, therefore, the assessed value of agricultural land. To 
the best of our knowledge, this is the first time that this 
trend has been demonstrated quantitatively with a predictive 
equation. 

The final set of predictive equations are those for change 
(delta) in assessed value 1967-72. Here it was found that 
changes in agricultural production do not explain changes in 


- nltUCUUMi’^' “ 
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Table 9 




ASSESSED VALUE 

AGRICULTURAL 

LAND* 



1967 



1972 


Indep.Var . 

B 

t ratio 

B 

t ratio 

Indep. Var 

V67001 

.490 

(5.21) 

.531 

(4.30) 

V72001 

V67008 

1.175 

(11.30) 

.910 

(7.78) 

V72008 

C 

3488.230 


6535.432 


C 


.707 


.562 


r2 

R2(a) 

.704 


.557 


R^(a) 

AVL1V67 

.632 

(5.89) 

.648 

(4.68) 

AVLIV72 

AVCRP67 

1.243 

(10.10) 

1.272 

(8.19) 

AVCRP72 

C 

3120.179 


5537.328 



r2 

. 695 


.589 



R^(a) 

.691 


.582 



AVLIV67 

,678 

(6.54) 

.701 

(5.35) 

AVLIV72 

AVCRP67 

1.181 

(9.89) 

1.179 

(7.94) 

AVCRP72 

SMSA67 

2986.172 

(3.01) 

5322.154 

(3.48) 

SMSA72 

C 

2191.150 


3340.624 



r2 

.725 


.639 



R^(a) 

.719 


.630 



Mean 

14545.002 


19894.229 



A Mean 


5349.227 





*A11 coefficients are significant at the. 01 level. 
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Table 10* 


TOTAL AGRICULTURAL ASSESSED VALUE 


1967 


indep.Var 

B 

t ratio 

V67001 

.949 

(6.08) 

V67008 

1.443 

(8.17) 

C 

8552.664 



.638 


R^(a) 

.634 


AVLIV67 

1.157 

(6.22) 

AVCRP67 

1.530 

(7.63) 

C 

7965.607 


r2 

.637 


R^(a) 

.633 


AVLIV67 

1.256 


AVCRP67 

1.396 


SMSA67 

6423.951 


C 

5967.050 


r2 

.699 


(a) 

.692 


Mean 

24880.412 


A Mean 


5681.38 


B 

1972 
t ratio 

Indep.Var 

.975 

(5.81) 

V72001 

.987 

(6.15) 

V72008 

12179.983 


C 

.546 


r2 

.541 


R^(a) 

1.164 

(6.24) 

AVLIV72 

1.404 

(6.71) 

AVCRP72 

10755.133 


C 

.580 


r2 

.575 


R2(a) 

1.234 

(6.97) 

AVLIV72 

1.281 

(6.39) 

AVCRP72 

7016.750 

(3.40) 

SMSA72 

7858.989 

- 

C 

.631 


r2 

.623 


R^(a) 


30561.794 


*A11 coefficients are significant at the. 001 level. 
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Table 11* 



DELTA AGRICULTURAL 

LAND VALUE (1967 

Indep. Var. 

B 

t ratio 

CHGLIV 

.607 

(2.87) 

CHGCRP 

.891 

(3.72) 

1021 

-1.357 

(8.07) 

1022 

.317 

(6.52) 

1023 

.089 

(5.70) 

SMSA72 

2330.842 

(3.08) 

C 

1817.123 



.561 


R2(a) 

.534 


Mean (10^) 

5349.227 



*All coefficients are significant at the .001 level. 
Except CHGLIV (.01) 
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Table 12* 


DELTA TOTAL AGRICULTURAL ASSESSED VALUE (1967-72) 


Indep. Var. 

B 

t ratio 

CHGLIV 

.792 

(3.05) 

CHGCRP 

.664 

(2.26) 

1021 

-2.072 

(10.03) 

1022 

.294 

(4.92) 

1023 

.165 

(8.61) 

SMSA72 

2299.777 

(2.47) 

C 

2603.019 


r2 

.639 


R2(a) 

.617 


Mean (10^) 

5681.383 



*A11 coefficients are significant at the .01 level except 

CHGCRP (.05) 

SMSA72 (.02) 
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assessed value. This is an expected result since we can assume 
that overall agricultural produccion in Ohio is probably near 
its peak and is only slightly affected by short term fluctua« 
tions in the produce market. We would then expect that it is 
the competition by other land uses for agricultural land which 
would better explain changes in assessed value. This is demon- 
strated by Tables 11 and 12. Here, one can see that changes 
in production combined with changes in the other land categories 
give a strong estimate of agricultural assessed value. 

There are actually several trends which seem to be measured 
by these equations. First, there is the effect of a slight 
increase in production. Second, there is the effect of urban 
land bidding up the prices of agricultural land. This is shown 
by the SMSA, commercial land, and residential land variables. 
Third, there is the effect of direct consumption of agricultural 
land on the urban fringe. This effect seems to be most highly 
correlated with the industrial land category and is why this 
coefficient is negative. Finally, it must be pointed out that 
changes in agricultural building assessed value are almost neg- 
ligible. The mean change in land value is $5349.227 x 10^ 
and in total value is $5681.383 x 10^ leaving only a change of 
$332.16 X 10^. This is because tax assessments between 1967 
and 1972 were made exclusively on the basis of market value of 
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land. The value of agricultural buildings for urban uses 
is negligible. Thus, the strongest equation for delta total 
agricultural assessed value is generated by using the same 
variables as those used for the land assessed value equation. 
This is shown in Table 12. 

Implications of the Agricultural Sector Models 

The agricultural sector models in combination with those 
for the other land sectors, allow estimation of the total tax 
base in Ohio counties. Given projections of population, 
employment, and agricultural production (input variables) the 
state will now have a tool to predict the resultant land con- 
version and tax base. 

Another important implication of these models is related 
to the usefulness of LANDSAT data in land use modeling over 
time. LANDSAT interpretations can more easily distinguish be- 
tween agricultural and urban uses than among urban uses. 

Since the urban and agricultural uses are shown to be inter- 
related by these models, it follows that LANDSAT can be used 
directly to monitor changes in agricultural resources due to 
urbanization and then, to help establish policies which will 
protect valuable agricultural land. The direct applications 
of LANDSAT for these purposes will be explored more fully in 
the final stage of this project. 
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raking all of the tax assessment and tax parcel models, 
alternative scenarios of future development in Ohio can be 
simulated and their land use tax base impacts can be estimated. 
In this way, public officials can better anticipate and plan 
for Ohio'p future. 
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Chapter II 


Converting Tax Base Data to Land Use Information 

One set of models which was produced during Phase I of this 
study is that for tax parcels. These models predict the number 
of tax parcels in each category of land use given projections 
of population and employment. Data on tax parcels are readily 
available and would, therefore, be an excellent, continuing 
source of information. What is needed then is an assessment 
of whether or not these data can be converted to actual land 
acreages. One of the major accomplishments during Phase II of 
the Ohio Land Allocation Model study is just such an assess- 
ment. The general approach, methodology, results, and con- 
clusions are explained in this chapter. 

Land Acreage f,.om Parcel Data 

In order to derive the acreage from data on parcels, one 
would need to know the size of each parcel. Alternatively, one 
may derive a frequency or probability distribution of parcel 
sizes in each county for each land category. This distribution 
can then be multiplied by the number of parcels in each cate- 
gory in order to get the acreage. Table 13 illustrates how one 
would calculate the probable acreage of residential land in a 
hypothetical county based on a sample of sizes of parcels in 
that county. 
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Table 13 

HYPOTHETICAL DISTRIBUTION OP RESIDENTIAL PARCEL SIZES 
SHOWING THE CALCULATION OP ACREAGE 


Parcel Size 
(Acre's ) 


Number in 
Sample 


Prequency Acreage * 


.1 

100 

.10 

250.0 

.25 

500 

.50 

3125.0 

.50 

200 

.20 

2500.0 

.75 

100 

.10 

1875.0 

1.00 

50 

.05 

1250.0 

1.50 

20 

.02 

750.0 

2.5 

10 

.01 

625.0 

7.5 

10 

.01 

1875.0 

IC.O 

5 

.005 

1250.0 

11.0 

5 

.005 

1375.0 

Total 

1000 

1.000 

14875.0 

*Actual total number 
size as example 

of residential 

parcels = 

25000 using 

25000 X .10 frequency = 2500 x 

.1 acres = 

250 acres 
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Sampling Parcel Data 


The purpose o£ the survey undertaken was to gain informa- 
tion about the size of tax parcels in each of the four major 
land uses emd to determine whether or not reliable profiles of 
parcel sizes could be established. This was done by tedcing a 
two percent random sample of the tcix parcels in three central 
Ohio counties: Pickaway CCircleville) , Delaware (Delaware), and 

Licking (Newark) . The land use type and acreage for each parcel 
selected were coded for keypunching and the samples were analyzed 
using the Statistical Package for the Social Sciences (SPSS) . 

The source of the parcel information was the public 
record of appraisal compiled by each county's auditor for 1975. 
Each peurcel record contained the owner's name, an identifying 
parcel number, a legal description of the p^.rcel, the assessed 
value of the land and buildings, and in some cases the acreage 
or dimensions of the parcel. These records are organized in the 
following way; 

Pickaway County 

First level - by township or incorporated area 
Second level - by school district 
Third level - alphabetically by owner's name 
Thirteen parcels per page 

Delaware County 

Same manner as Pickaway, but 35-40 records per page 


Llclclng County (two complete listings) 

Alphedsetically by owner's name, thirteen records 
per page 

First level - by township or incorporated area 

Second level by school district 

Third level - by land use type 

Fourth level - alphabetically by owner's name 

Sixteen records per page 

A random number table with values ranging from 1 to 52 was 
used to select the parcels to be recorded. In Pickaway and 
Licking Counties the rauidom number between 1 and 52 identified 
one record in a four page block. In Delaware County, only values 
between 1 and 40 were used and each number identified a record 
for a single page. When a parcel was selected, the land use and 
acreage of the parcel was recorded with the following exceptions: 

1) If the parcel had no acreage recorded, its sub- 
division name euid lot number were noted and the 
dimensions of the lot taken from the plat record, 
converted to acres and recorded. 

2) If the legal description contained the dimensions 
of the parcel, e.g. 150' x 250', this information 
was noted and converted to acres. 

3) If the legal description did not identify a sub- 
divison name or contain the dimensions of the 
parcel, the parcel was rejected and another random 
number and parcel selected. A deed search is 
necessary to identify the acreage of such a parcel. 
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Because industrial and commercial parcels make up a very 
small portion of the parcels within a county# very few such 
parcels were selected in the random sampling process. This de- 
ficiency was corrected by selecting a township book in an area 
containing industrial and commercial land and recording all such 
parcels which showed acreage in the sample until a minimum of 
fifty parcels were recorded. 

The coded information was keypunched and SPSS was used to 
calculate the frequency distributions and simple statistics for 
each sub-sample# in both grouped and ungrouped form. Table 14 
displays the results of the ungrouped analysis and the sample 
sizes for each land use and county. 

The frequency distributions obtained from the random 
sample were compared using the Smirnov test. This test was 
chosen because it does not assume a normally distributed popula- 
tion# and the populations from which the samples were drawn are 
not normally distributed. The test statistic, called the D value# 
is the largest difference between the cumulative frequencies for 
each grouping of the two samples being tested. The samples were 
grouped as shown in Table 15. If the D value is larger than the 
critical value, the two samples are significantly different. The 
samples were compared within land use types, the results are dis- 
played in Table 16. The significantly different samples were the 
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Table 14 

PARCEL SIZE SURVEY SIMPLE STATISTICS {UNGROUPED DATA) 


All 


Agricultural Parcels 

Pickaway 

Delaware 

Licking 

Counties 

N 

106 

156 

235 

497 

Mean 

49.3 

40.4 

43.4 

43.7 

Median 

32.0 

38.1 

32.8 

34.0 

Variance 

2791.3 

1038.8 

1523.9 

1644.6 

Skewness 

1.449 

1.668 

1.252 

1.557 


Industrial Parcels 





N 

53 

73 

81 

207 

Mean 

11.3 

9.0 

23.6 

15.3 

Median 

3.0 

3.1 

6.0 

4.2 

Variance 

413.0 

173.3 

1560.1 

163.5 

Skewness 

3.394 

2.038 

2.285 

6.929 

Commercial Parcels 



* 


N 

52 

80 

176 

308 

Mean 

3.0 

3.5 

4.8 

4.2 

Median 

0.8 

0.8 

0.9 

0.8 

Variance 

139.3 

49.9 

222.6 

163.5 

Skewness 

6.628 

3.429 

6.457 

6.929 


Residential Parcels 
N 

240 

530 

754 

1524 

Mean 

1.0 

1.4 

0.8 

1.0 

Median 

0.3 

0.3 

0.2 

0.25 

Variance 

4.2 

8.5 

3.6 

5.5 

Skewness 

4.404 

3.643 

6.514 

4.783 


All Parcels 




N 

451 

339 

1246 2536 

Estimated Total 

Parcels 18,000 

29,000 

52,000 99,000 

Sample Size 

2.5% 

2.9% 

2.4% 2.6% 

♦State Board of 

Tax Appeals 


, is 
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TABLE 15 




Land 

Use Type 


Group 

Agricultural 

Industrial 

Commercial 

Residential 

1 

.1-.9 

.1-.9 

.l-;2 

,1 

2 

1.0-9. 9 

1.0-4. 9 

.3-. 4 

.2 

3 

10.0-29.9 

5. 0-9. 9 

. 5-.6 

.3 

4 

30.0-49.9 

10.0-14.9 

.7-. 8 

.4 

5 

50.0-69.9 

15.0-24.9 

.9-1.0 

.5 

6 

70.0-89.9 

25.0-49.9 

1.0-2. 9 

.6 

7 

90.0-109.9 

50.0-74.9 

3. 0-4. 9 

.7 

8 

110.0-129.9 

75.0-99.9 

5. 0-6. 9 

.8 

9 

130.0-159.9 

100.0-149.9 7. 0-9. 9 

.9 

10 

160.0+ 

150.0+ 

10.0-49.9 

1.0-1. 4 

11 



50.0-99.9 

* 

1 

12 



100 . 0+ 

2. 0-2. 9 

13 




3. 0-4. 9 

14 




5.0+ 


Inclusive Group Limits in Acres 
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Table 16 

SUMMARY OF SMIRMOV TESTS 




Samples Tested 

AGRICULTURAL 

LAND USE 
INDUSTRIAL 

TYPES 

COMMERCIAL 

RESIDENTIAL 

Licking-Delaware 

N{1) 

N(2) 

N(l) N(2) 

N(l) N(2) 

N(l) 

N(2) 

• 

235 

156 

81 73 

176 80 

754 

530 

Critical Values 







.05 Level 


.140 

.219 

.184 


.077 

«01 Level 


.168 

.263 

.220 


.092 

D Values 


.091 

.184 

.216* 


.158* 


Pickaway-Delaware 

N{1) N(2) 

N(l) N(2) 

N(l) N{2) 

N(l) 

N(2) 

Critical Values 

106 156 

33 73 

52 80 

240 

530 

.05 Level 

.171 

.245 

.242 


.106 

.01 Level 

.205 

.294 

.290 


.127 

D values 

.139 

.184 

.135 


.094 


Licking-Pic)'away 

N(l) N(2) 

N(l) N{2) 

N(l) N{2) 

N(l) 

N(2) 

Critical 'lues 

235 106 

81 53 

176 52 

754 

240 

.05 L 1 

.159 

.240 

.215 


.101 

.01 Level 

.191 

.288 

.257 


.121 

D Values 

.102 

.061 

.081 


.215* 


♦significantly different samples 


Critical Value at 
Critical Value at 


.05 Level 
.01 Level 


1 . 3 ^ 

1 . 63 '\ 


N(l) + N(2) 
'n(I) X N(2) 
HTd) + NviT 


,N(1) X N(2) 


D Values = Maximum (Cumulative Frequency (1, i) - Cumulative Frequency (2, i) ) 
i = 1 to number of groups in distribution 
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residential sample for Licking County compared to both Delaware 
and Pickaway Counties, and the commercial samples for Delaware 
and Licking Counties, The difference probably resulted from the 
more urbanized nature of Licking County. Platted parcels are 
more numerous in urbanized areas and are generally smaller than 
non-platted parcels. 

Conclusions 

The three samples taken do not allow generalization of 
the results to the entire state. However, they do indicate that 
profiles of different types of counties (urban, rural, metro- 
politan) could be developed with a sample of 10-15 counties. 

Sampling this many counties could require as many as 1000 man-hours. 
However, as more auditors* records become computerized, such a 
sample could become feasible. If reliable profiles could be de- 
veloped it would allow translation of the existing number of 
parcels data to acreage, providing much useful information about 
land use and land use change. In this way, projections made 
using the tax parcel models can be directly translated into land 
use projections. 

The results of this type of analysis could then be used as a 
continuing and consistent check on the accuracy of macro level land 

I 

use totals produced by LANDSAT. Future work with the tax models and 
continued use of LANDSAT could generate a data base in the future which 
directly linked land use change as measured by LANDSAT with the changes 
in tax base in Ohio Counties. 

niVlonUCtBILiTY OF 1 * 
li.hriNAL PAGE IS poor. 
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Chapter III 


Land Conversion in Franklin County 
Establishing a Test for LANDSAT Data 

LANDSAT offers the opportunity to provide computer com- 
patible land use information at frequent intervals. Thus, land 
use changes as measured by LANDSAT can be linked with changes 
in population, employment, and other socio-economic character- 
istics to produce empirical models of land use cheuige. This is 
the underlying goal of the modeling effort in the State of Ohio. 

In pursuing this goal, a number of technical problems have 
arisen which require further exploration before a set of final 
models can be produced. This chapter summarizes the work using 
aerial photography for Franklin County, Ohio, which will be 
utilized in performing a test of the accuracy of LANDSAT inter- 
pretations . 

Land Use Modeling and LANDSAT 

In order to create a predictive land use model for counties 
in Ohio , a land use data base is required which is both extensive 
and accurate. LANDSAT provides a unique opportunity for such 
modeling because it is the most extensive, most consistent data 
base available. Unfortunately, a number of technical problems 
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arise in converting from I>ANDSAT imagery to land use categories. 
These problems leave a potential for large errors in the final 
land use classifications. Such errors would be unacceptable 
for the purpose of land use modeling, i.e. relating land use 
changes to socio-economic cht.nges. Thus, the first step in 
analyzing LANDSAT must be to identify the nature and degree of 
these errors, and if possible, to derive correction factors 
which might be applied prior to the use of these data in a 
model . 

The potential errors associated with LANDSAT have been 
subdivided into two major components. Those errors associated 
with the misalignment of pixels for two different LANDSAT scenes 
are referred to as Error 1. This error results from the pos- 
sibility of + one pixel misalignment in ground orientation as a 
result of the rescanning readjustment, and reclassification of 
the original, distorted data pixels. 11 The second type of 
error (Error 2) results from the misclassification of land 'cover 
due to the similarities in spectral signatures of different land 
uses. To date, the nature and extent of these errors have not 
been quantified. This is the first task of LANDSAT data analysis 
in Ohio. 

^^Space Applications Board, Assembly of Engineering, Nat. Res. 
Cotoncil, Supporting Paper 3 Land Use Planning (Washington, 

D.C. Nat. Acad. Sci., 1974), pp. 31-32, 

^^Bendix Aerospace Division, Computer Mapping of LANDSAT Data 
for Environmental Applications (Ann Arbor, Michigan: Bendix 

Corporation , November 1975 )' , p . 8 . 
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Quantifying LANDSAT Errors 

In order to quantify the misclassif ication errors as- 
sociated with LANDSAT, a number of tests have been devised. 

The first test relates to the overall changes in land use 
which can be tabulated from LANDSAT data. Scenes for Franklin 
County for both 1973 and 1975 will be interpreted by Bendix 
Corporation and given to the State of Ohio on computer tapes. 
Based on these tapes# one can calculate the changes in land 
use in Franklin County during this period. The question which 
arises is how accurate are these land use change figures? 

Since LANDSAT has been in use for such a short time# no one 
appears to have utilized the data in this way or to test the 
accuracy of the results. 

Fortunately# two sets of aerial photographs are avail- 
able for Franklin County at times very close to those for 
LANDSAT imagery. Thus# the first step in compiling information 
for a test of LANDSAT involved the interpretation of land use 
change from these aerial photographs. 

The initial task was to correlate the 1972 Franklin County 
aerial photographs (scale 1:1000) with the 1976 set (scale 
1:2000). Upon completion and matching of the comparable mosaics 
with differing scales# the photographs were scanned on a zoom 
transfer scope for land use changes . Any differentials apparent 


reproducibility or I 
'ORIGINAL PAGE IS PCv. 

- 37 - 


in the two sets of photos were noted on the 1972 photos. We 
attempted to identify the entire area of change during this 
process , noting not merely the existence of a new structure 
but also its attendant property lines, if they could be iden- 
tified from the landscaping. As we were interested in land 
use changes only, differing crop patterns and water levels 
were considered extraneous. 

Categorization of land use changes appeared under seven 
classifications; 1) Urban Recreation 2) Open Space - all land 
void of development that could not be identified as agri- 
cultural land 3) Agriculture 4) Commercial - including trucking, 
airports, and warehousing operations 5) Industrial 6) Resi- 
dential 7) Public - including schools, churches, and highways, 

A separate notation of "A" was utilized to signify land cleared 
and/or under construction and would be appended to the proper 
land use symbol. 

Upon comparison of the areas on the 1972 set of aerial 
photos with the actual changes in the 1976 photos, each land use 
change was catalogued as outlined above, with the 1972 land use 
recorded a separated from the recorded 1976 land use change by a 
slash (e.g. 3/6). Certain decision rules were derived for clas- 
sification of land uses as follows; 
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1) Commercial buildings - a higher parking lot to building 
ratio than industrial structures. Location along a 
major arterial. Location as related to residential de- 
velopment. 

2) Industrial - location near interstate interchanges . 

Parking lot to building ratio lower than commercial. 
Landscaping and large front set back. 

3) Additional buildings on a previously established land 
use was not noted, except in the case of new resi- 
dential units in subdivision development. 

4) Individual parcels in subdivisions cleared for deve- 
lopment in 1972 and developed by 1976 were noted as 
6A/6, 

5) Agriculture - clear crop pattern with evidence of 
farm buildings in close proximity. 

Once the .classification was completed and noted on the 1972 

photos, a planimeter was j.ied to determine the areal change. 

Thus, a table was derived detailing for each 1972 photo the land 

use changes on the photo with its before/after classifications, 

the area of land use change in square inches, and finally, this 

square inch measurement converted into acres. Table 17 shows 

these data. 

In order to quantify the errors associated with LANDSAT 
at the pixel level, we made use again of the aerial photographs. 
First, the photos were oriented to the U.S.G.S. quadrangles using 
a zoom transfer scope. Then, a 1.1 acre grid of the same dimen- 
sions as LANDSAT pixels was overlayed on the mosaic of photos for 
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Table 17 

LAND USE ClfANGES IN FRANKLIN COUNTY BETWEEN 1972 AND 1976 

(In Acres) 


Urban Open A* 

Recreation Space Agriculture Coitangrclal Industrial Residential Residential Public Total 


A* 

Urban 

Recreation 

Open 

Space 


4.591 


136.823 


Agriculture 344.582 
A* 

Agriculture 


Conmercial 


A* 

Commercial 


Industrial 


A* 

Industrial 


Residential 


A* 

Residential 


A* 

Public 


485.996 


16.750 


5.510 


8.724 


1^154.270 


366.850 


196.281 


68.871 


€37.052 4,210.055 

709.826 2^278.466 


220.156 


14.233 


2^542.470 


14.234] 


1.786.272 


1,581.267 


9,030.991 


328.053 


219.927 


547. 9B0 


6B6.869 


61.983 


4.591 


20.202 


367.080 


1,140.725 


3 , 981.634 


4.591 


196.201 


5.510 


220.156 


112.030 


367.080 


14 , 604.223 


A* - under construction and/or land cleared for development 

Sourcei 1972 fFlight Ho- 5010) and 1976 (Flight Ho. 5979) aerial photos of Frankllii County with field checks. 

























i 

I 

I 

t 

each year, for a sample quadrangle in Franklin County. For 
each "pixel" the dominant land cover (greater than 50%) and 
the probable land use, if different, was recorded, encoded, 
and punched onto computer cards* In this way, a set of 
highly accurate land use information was produced which is 
compatible with the format of the LANDSAT data. When the 
LANDSAT data become available, a computer generated pixel 
by pixel comparison may be made of land cover differences, 
and a quantified index of Error 2 can be generated. 
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Chapter IV 

Study Design for the Final Phase 

During the final phase of the project, a number of addi- 
tional tasks will be accomplish^. The first of these is the 
tests of the accuracy of LANDSAT data. Error 2 will be quan- 
tified by comparison with the aerial photograph "pixel” data 
described in Chapter III. 

Identifying Error 1 is a much more difficult task. Ideally, 
one would take two scenes of the same area with a small time gap 
between them (e.g, 18 days) . Any changes in land use recorded 
for this time gap would have to be due to misalignment since a 
real change is unlikely. Unfortunately, different weather con- 
ditions, sun angles, etc. can produce two distinct sets of 
spectral signatures in eighteen days. Thus, errors identified 
over this period would be due not only to misalignment but also 
misclassification and changes in the training set. For the pur- 
poses of this study then, we will analyze a one day "gap." This 
will consist of the one day side lap area for Franklin County. 

In one day, spectral signatures are not likely to change. Yet 
problems with aligning these two different scenes will occur 
and should be reflected in this overlap area. To the extent 
that this overlap is not a true measure of the alignment problems 
which occur from year to year, we will not identify the mis- 
alignment error. However, it is felt that this analysis will 
shed some light on this problem. 
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Correcting LAMDSAT Errors 


Given a determination of the errors associated with LANDSAT, 
the next step will be to see if these errors are consistent. For 
example/ it may be that certain types of agricultural land are 
consistently identified as urban. There may be a definitive per- 
centage of error in urban categories across all types of errors. 
Given this information, error correction factors can be derived 
and applied to all LANDSAT data in Ohio. 

Not only can the aerial photo data be used to develop cor- 
rection factors for the LANDSAT data provided by Bendix Corporation, 
but they can also be utilized to test the relative accuracy of other 
methods of LANDSAT computer compatible tape data analysis. Thus, 
one could reinterpret the Franklin County LANDSAT tapes with addi- 
tional, new, or revised techniques and algorithms utilizing the data 
gathered in the course of this project as the baseline, "real world" 
measure of accuracy. 

It would be expected that all methods of interpretation would 
base some level of inaccuracy for the reasons discussed above. How- 
ever, the method of analysis and data presented in this report would 
help to delineate the level of error associated with a number of 
alternative interpretation techniques. 

Deriving a Land Use Model 

Once error correction factors have been derived and applied 
to LANDSAT, the data can be utilized to generate a set of land 
use models. Regression analysis will be used to determine the 
amount of land in each county in each use as a function of socio- 
economic characteristics in the same time period. The nature of 
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these models will be the same as those derived in Phases I and 
II of this project. 


i 
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Conclusions and Model Outputs 

Having performed all of the tasks briefly outlined above, 
a number of useful outputs will be produced. First, a quan- 
titative definition of the errors associated with current 
LANDSAT interpretation will be made. This can be used to point 
to areas where further technical development is required. 
Second, the State of Ohio will have a measure of the accuracy 
of its LANDSAT interpretations and a set of error correction 
factors. Finally, LANDSAT data will be linked in a modeling 
framework with social and economic data. This model can be 
utilized in land use planning efforts at the state, regional, 
county, and local levels. 

It will then remain to supervise the programming and 
mounting of each of the computer models produced during the 
cotirse of the study, on the state computer system. Table 18 
shows models on tax base produced in Phase I of the project. 

In addition to these 34 models, there will be three additional 
models of agricultural assessed value and at least one model 
of LANDSAT land use change (assiiming the accuracy level is 
acceptable . 

Each of these models will be documented as follows; 

1. Explanation of what the model predicts. 

2 . Explanation of the required input variables . 

3. Specification of formats for input variables in card deck. 
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Table 18 


PROPOSED GENERAL FORM FOR T.AX MODELS 


Potential Models 


Static Models 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


Change Models 


Residential Parcels 18 

Commercial Parcels 19 

Industrial Parcels 20 

Agricultural Parcels 21 

Residential and Commercial Parcels 22 ' 
Residential AV Land 23 

Residentia'. AV Total 24 

Commercial AV Land 25 

Commercial AV Total 26 

Industr al AV Land 27 

Industrial AV Total 28 

Agricultural AV Land 29 

Agricultural AV Total 30 

Residential and Commercial AV Land 31 


Residential and Commercial AV Total 32 
All Taxable Land AV 33 
All Taxable Total AV 


34. 


For a run of one or more models, the user would prepare a card 
deck containing the following data: 

Card 1 - number of counties to be processed 
Card 2 - name of 1st county 

Card 3 - munber of models desired (NEQS) and the index 
number of each model equation 

Card 4 through (3 + NEQS) - input values for each 
model requested 

End of file marker or repeat of Cards 2 through end 
for each additional county. 

In this way, one can forecast the tax base, parcel numbers, 

* 

and land use for up to 88 counties in Ohio utilizing input pro- 
jections of population, employment, agricultural production, 
and the SMSA dummy variables. One can select to use only those 
models which are relevant to a particular project. Alternative 
projections of population, employment, and the other input 
variables can be used to simulate the potential impacts on the 
property tax base and land use. 

These models can also serve as control totals for land use 
and tax changes at the county level. Then, further model deve- 
lopment may be carried out at the micro level to yield additional 
allocation procedures. All of the models should thus serve to 
aid the work of public officials and planners at the local, 
regional, and state levels. 
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APPENDIX A 


AN ANALYSIS OF THE PEBFOPMANCE OF THE DEMOS MODEL 
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Introduction 

The major purpose of the State's land use modeling effort is 
to find a functional relationship between variables to land use 
and land use change (dependent variables) and variables related 
to economic and demographic characteristics of the state (ex- 
planatory variables independent variables) . The explanatory 
variadiles Include employment by category, population size, and 
the changes in those variables over time, at the county level. 

The report on Phase I details the statistical analysis which was 
required to define the mathematical relationships among the 
dep.endent and independent variables. 

Given a forecast of the economic and demographic variables 
and using these mathematical equations, it is possible to predict 
the new patterns of land use, land use changes, and thus, to 
analyze the level of pressure of land use conversion in each 
county. The accuracy of the prediction will depend not only on 
the specification and statistical performance of the functional 
equations in the land use model, but also on the performance of the 
State's model for forecasting economic and demographic variables 
at the county level. It is for this reason that we were asked to 
undertake an analysis of DEMOS, the model the State currently uses 
for these projections. The results of this analysis are summarized 
below. 
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The DEMOS Model 

DEMOS# a model developed by Battelle Memorial Institute, is 
the model used by the State to forecast economic and demographic 
variables. Given two key input parameters the national unemploy- 
ment rate and a population growth rate (assumed to be either slight 
or no growth) ,• DEMOS predicts employment in thirty-nine categories 
and population by age group, all at the county level. The projec- 
tions are ntade for each county taken one at a time and can be 
derived for all 88 f^Mo counties. The starting date or base period 
for the model is 1970 and the projections can be made on an annual 
basis to 1985. Table A-1 shows an example set of projections for 
one Ohio county. Because population change variables are not used 
in the land use change model equations, only the DEMOS employment 
projections were analyzed. 

Testing Performance of DEMOS 

For the four year period 1970-1974, DEMOS projections can be 
considered to be estimates of employment in each county. Thus, 
estimates were made for all 88 counties in Ohio for 1972 and 1974 
using actual national unemployment rates for these years and as- 
suming the zero population growth rate which has occurred over the 
past few years. These estimates can be compared with real counts 
of employment by county dtiring the same period available from the 
Ohio Bureau of Employment Services (OBES) . An analysis of the 
performance of DEMOS can be made if appropriate aggregation of 
DEMOS categories into OBES categories is carried out. 
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TABLE A-1, Example DEMOS Forecast 


POPULATION BY 5- YEAR AGE GROUP, RICHLAND COUNTY, 1972 


TOTAL 

133,596 

0-4 

11,659 

5-9 

12,700 

10-14 

13,959 

15-19 

12,885 

20-24 

10,696 

25-29 

9,447 

30-34 

8,455 

35-39 

7,737 

40-44 

7,780 

45-49 

7,787 

50-54 

7,183 

55-59 

6,166 

60-64 

5,170 

65-69 

4,156 

70-74 

3,162 

75 + 

4,654 




EMPLOYMENT BY INDUSTRY, RICHLAND COUNTY, 1972 


TOTAL 51,926 
AGR, FORESTRY & FISHING 906 
MINING 55 
CONSTRUCTION 2,306 
FURNITURE, LUMBER, WOOD 253 
METALS INDUSTRY 5,650 
MACHINERY, EXCEPT ELECTRICAL 1,796 
ELECTRICAL MACHINERY 6,857 
TRANSPORTATION EQUIPMENT 1,497 
OTHER DURABLE GOODS 2,860 
FOOD AND KINDRED PRODUCTS 216 
TEXTILES AND TEXTILE PRODUCTS 197 
PRINTING AND PUBLISHING 1,104 
CHEMICALS 80 
OTHER NONDURABLE GOODS 2,381 
RAILROAD AND RAILWAY EXPRESS 172 
TRUCKING 979 
OTHER TRANSPORTATION ISO 
COMMUNICATIONS 857 
UTILITIES AND SANITARY SERVICE 437 
WHOLESALE TRADE 1,690 
FOOD AND DAIRY STORES 1,098 
EATING AND DRINKING PLACES 1,714 
GENERAL MERCHANDISING 1,802 
MOTOR VEHICLE RETAILING 965 
OTHER RETAIL TRADE 2,089 
FINANCE 563 
INSURANCE AND REAL ESTATE 1,502 
BUSINESS AND REPAIR SERVICES 1,236 
PRIVATE HOUSEHOLDS 491 
OTHER PERSONAL SERVICES 1,429 
ENTERTAINMENT 301 
HOSPITALS 1,335 
OTHER HEALTH SERVICES 920 
GOVERNMENT EDUCATION 2 ,'260 
PRIVATE EDUCATION 528 
OTHER EDUCATIONAL SERVICES 176 
RELIGIOUS AND NONPROFIT ORGS 653 
PROFESSIONAL ORGANIZATIONS 675 
PUBLIC ADMINISTRATION 1,716 
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The DEMOS output was aggregated down to seven categories 

equivalent to those of OBES. These are shown in Table A-2 being 

in general terns four basic activities and three non-basic. 

TABLE A- 2. Bureau of Employment Services (OBES) 

Employment Categories 

Mining Wholesale and Retail Trade 

Contract Construction Finance, insurance and 

Manufacturing Real Estate Services 

Transportation and 
Utilities 

Comparisons were made for two sub-periods, 1970-1972 and 1972- 
1974 in order to test the sensitivity of the DEMOS model to 
national trends. 

It was recognized from the outset that DEMOS projections of 
employment are based on place of residence while those of OBES 
are based on place of work. This made direct comparison of 
category by category employment figures for each county impos- 
sible. However, the changes in each set of data (e.g,, OBES 
1970-72; DEMOS 1970-72) should exhibit exactly the same trends 
if the DEMOS employment projections are accurate. 

The statistical analysis performed was to regress in a 
cross-sectional way, the change (A) employment for the sub- 
period (1970-72 or 1977-74) , to employment at the start of the 
sub-period. Functionally, it is expressed as; 
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(70-72) = a| 

^i 

(70) 

r2 

^d 

<I 

(72-74) = A? 

JU 

yd 

^i 

(72) 

r2 

^d 


where 

= change in (Delta) employment in category i 
given by DEMOS 

= employment in category i given by DEMOS at 
the start of the sub-period 

i = employment categories (1, 2, 3,...7) 

and are defined below. 

# of observations = 88 

Two statistical parameters are important in this analysis: 

1) The regression coefficient A^^ can be interpreted as 
the state’s average growth rate of the economic 
sector i or in other words, how much of the growth 
is common to all the counties independent of their 
internal and local attributes. By internal and 
local attributes we mean such things as industrial 
attractiveness, agglomeration economies and economies 
of scale, all of which vary widely from county to 
cotinty . 

0 2 

2) The coefficient of determination (Rjj) can be inter- 
preted as how much of that average growth in the 
sector is tied to the state's growth. The lower the 
r 2, the more important are the local attributes in 
explaining growth and thus changes in population and 
employment. The larger the r 2, che less important 
are these local attributes. In general, we can ex- 
pect that for the basic sectors, the local attributes 
are important, i.e. we will obtain a lower Br . 
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Results of Analysis 


SECTOR 


1970 - 

1972 

1972 - 

1974 

1 


OBES 

DEMOS 

OBES 

DEMOS 


A 1 

.216 

.008 

-.009 

1 

o 

o 

00 

Mining 1 , 



.54 

.92 

.002 

.75 


A i 

-.007 

.004 

-.0006 

-.017 

Const. 2 . 


R^ 

.002 

.92 

. 00004 

,87 


A 1 

-.086 

.027 

.029 

.026 

Manuf . 3 „ 


R^ 

.84 

.96 

.40 

.996 


A 1 . 

-.035 

.004 

.004 

.015 

Transp . 4 


.39 

.80 

.02 

.88 


A i 

.067 

.007 

.060 

.014 

Wholesale 5 

r2 

.71 

.34 

.85 

.71 


A i 

.102 

.007 

.039 

.014 

Financial 6 

r2 

.81 

.66 

.66 

.75 


A i 

.773 

.007 

.097 

.015 

Services 7 ^ 


R‘ 

.99 

.51 

.98 

.75 


A i 

.066 

.013. 

.048 

.018 

Total 

r2 

.70 

.87 

.88 

.91 


The same analysis was performed on the OBES data. 


Xj C70-72) = a2 X° (70) r2, 

X9 (72-74) » A? X? (72) r2 

where the symbols represent the same parameters for the OBES data. 

The A®'s and R^'s of this analysis can be used as benchmark 
figures to compare against the Aj^'s and of the previous step. 

In this way the performance of DEMOS can be analyzed. 

The results of the statistical analysis show a large dis- 
crepancy between DEMOS estimates and OBES benchmark figures. This 
is demonstrated in Tables A-3 and A-4. As an example we can take 
the manufacturing sector. In the sub-period 1970-72, OBES shows 
a substantive decrease (A = .086; A^ = ,029) in manufacturing 

employment against a gain by DEMOS . Comparing the R^ values , we 

2 

see that OBES shows that on 40% (R^ =40) of the variation of 
the growth rate all of the growth of the sector (R^ = 996 or 
99.6%) is common across the counties without any consideration of 
local conditions. 

Conclusions 

Comparison of DEMOS projections for the period 1970-74 with 
actual employment data for the same period has shown the model 
to make substantial errors. Population projections could not be 
analyzed because no actual population counts are available for 
the years other than 1970. Use of these projections to simulate 
land use changes based on the equations derived in the present 
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project would thus lead to erroneous expectations with regard 
to land use change. Thus, the use of the DEMOS model employ- 
ment projections for this purpose is not recommended. 
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TABLE A-4 


Differences Between DEMOS and OBES Predictions 

1970 - 1972 1972 - 1974 



OBES 



DEMOS 

OBES 


DEMOS 


Mining 

.216 

.54 

-.008 

.92 

.009 

.002 

-.008 

.75 

Constiuctlon 

o 

0 
• 

1 

.002 

.004 

.92 

-.00006 

.00004 

.017 

.87 

Manufacturing 

-.086 

.84 

.027 

.96 

.029 

.40 

.026 

.996 

Transportation 

-.035 

.39 

.004 

.80 

.004 

.02 

.015 

.88 

Wholesale 

.067 

.71 

.007 

.34 

.060 

.85 

.014 

,71 

Financial 

.102 

.81 

.007 

.66 

.039 

.66 

.014 

.75 

Services 

.773 

.99 

.007 

.51 

.097 

.98 

.015 

.75 

Total 

.006 

.70 

.013 

.87 

.048 

.88 

.018 

.91 


G R 

r2 

G R 

r2 

G R 

R^ 

G R 

r2 
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APPENDIX B 

The Newlinq Model: Testing a Method For 

Forecasting Population for Minor Civil Divisions 

By 

Harvey Curran 
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INTRODUCTION 


The purpose of this study was to find a quick and dirty method 
of forecasting future population for minor civil divisions (MCD's). 
Population forecastina seems important because it forms the basis 
for much long-range planning. Quick and dirty methods are desirable 
because the time and expensive computing machinery to do rigorous 
work may not be available to a planner who needs a forecast. Minor 
civil division projections are desirable because available forecasts 
are for SMSA's or counties rather than MCD's. A survey of population 
forecasting literature led to the method developed at Rutgers University 
by Michael Greenberg and others using the Newling density model. ^ 

The Greenberg-Newl ing method is appealing for several reasons. 

It requires little input data, the calculations can be performed on 
any calculator capable of exponentiation, and the method allows the 
forecaster to incorporate his/her knowledge of special conditions. The 
major question about the method is its accuracy. To test the method, 

I followed the Greenberg example for calculating the model parameters 
using Ohio data and applied the results to each townshio and community 
in the Miami Valley planning region. This region was chosen because it 
includes a representative variety of rural, suburban, and urban areas. 
The following sections present the Newling model, the procedures used to 
calculate the parameters, and the results of applying the model to the 
Miami Valley region. 
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THE NEWLIN6 MODEL 

This entire section is taken directly from Seciton 4.2 of the 
Greenberg book. 2 

The Newling method involves the derivation of critical population 
densities for urban, suburban, and rural MCD's. All communities are class- 
ified and will converge to one of the critical densities, some like older 
central cities by losing copulation, others like former farmlands adjacent 
to central cities by gaining population. Eventually all of the MCD's in 
the region will approach a steady state of clearly defined rural, suburban, 
or urban densities. 

Newling originally developed his model after observing an inverse 
relationship between relative population growth from 1950 to 1960 and 
pooulation density in 1950. This relationship is expressed in Equation 2.1, 

Eq. (2.1) (l+r^^) = Ad^k 

where (l+»'^ ) is the growth ratio, i.e. the rate of 

t growth at each density for a given density 
ceiling group, 

A is a constant for each density ceiling 
group. It is the growth ratio when the 
density is one person per unit area, 

K is a constant for each density ceiling 
group. It is the ratio of the rate 
of change of growth to the rate of change 
of density, and 

dt is the density at time t. 

The density at one time period forward (d^+-|) may be expressed as 
Eq. (2.2) . (1 t , 

or, = Ad.j.“k(j^ 

which may be simplified as 
Eq.. (2.3) dt+l = Adtl-k . 

Equation (2.3) may be generalized for time (t+m) where m is any number 




I 


I 

i 

r of intervals. 


(1-k)'" 


Eq. {J?.4) dt+m 


^1/k (^1/kj n-k)"’ 


The expression within the large brackets approaches unity as (m) 
becomes very large. Thus is the critical density. Each MCD's 

future densities are related to its initial density (dt). the critical 
density of its density ceiling group, and the rate (k) at which 

all members of the group approach the critical density. 

Newling determined the parameters (A) and (k) for three groups 
of counties by least squares. All communities were classified with 
reference to their density in 1950 and growth rate from 1950 to 1960 by 
partitioning the scattergram into segments equally removed from the 
regression lines. The density of each unit is projected by substitutina 
into Equation (2.4) the appropriate (A) and (k) valuer and its pooulation 
density in 1960 as the value of (d^). 

In essence, density acts as a surrogate for a host of previously 
cited factors responsible for the suburbanization process. For example, 
the procedure can subsume the decline of old densely developed central 
cities, the rapid increase and then the leveling off of suburbs, the 
preservation of lands in rural land uses, and other conmonly observed 
phenomena. 


DEVELOPING (A) AND (K) PARAMETERS FROM OHIO DATA 

Following the example of Newling and Greenberg, We used least squares 
regression to determine the (A) and (k) parameters for rural, suburban, 
and urban density ceiling groups using data from Ohio counties. The first 
step in the orocess is to make a scattergram of the natural log of the 1950 
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density for each Ohio county against the natural log of the average decade 

growth rate between 1950 and 1970 for each county. The average decade 

growth rate is defined as 

Eq. (3.1) GROWTH * 1970 Population / 1950 Population 

If GROWTH > 1 then GROWTH = (GROWTH-1 )/2 + land 

If GROWTH < 1 then GROWTH = (l-GR0WTH)/2 + GROWTH. 

The resulting scattergram was examined for linear groupings and 
divided into three groups. Least squares regression was performed on 
these three data groups with outlying points being eliminated to improve 
the fit. The regression equations which resulted and their R^ and F 
values are displayed in table 3.1. 


TABLE 3.1 Regression Equations Developed from Ohio Data 
Log (GROWTH) = a + B*Log(1950 Density) 


Class 

a 

b 

r2 

F 

Rural 

.97718 

-.22889 

.71 

37.8 

Suburban 

1.19257 

-.24196 

.86 

84.6 

Urban 

.97767 

-.1167 

.90 

102. f 


The parameters (A), (k), and (A^^**^) were developed from the regression 
coefficients using the relations in Equation (3.2). 

Eq. (3.2) (A) = antilog(a) 

(k) = -B. 


R15PROt)UCTrV* 
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Table 3.2 displays the (A) and (k) parameters for each density 
celling class along with the critical densities. 

TABLE 3.2 (A) and (k) Parameters .and Critical Density 


Class 

A 

k 

Critical Density 

Rural 

2.657 

.22889 

71.47 (persons/sq. mi.) 

Suburban 

3.2955 

.24196 

138.21 

Urban 

2.6583 

.1167 

4348. 70 


The regression equations were graphed and the graph divided into 
sections for use in classifying the MCO's to be forecasted. Figure 
3.1 displays the subdivided graph of Log(1950 Density) against Log(GROWTH). 
The solid lines are the regression lines and the dotted lines are the 
dividers between density ceiling classes. 
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FIGURE 3.1 Subdivided Graph for Determination of Density Ceiling Classes 
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TESTING THE NEWLING METHOD 


The inout data necessary for the Newling model was co^ected from 
the 1950 and 1970 Censuses of Population and the Ohio Department of 
Natural Resources OCAP system. The input data set consists of the 1950 
population, the 1970 population, and the area of each MCD. The first 
step in applying the model is the classification of MCD's by density ceiling 
classification. The 1950 density in persons per square mile and the 
average decade growth rate were calculated, the natural log of each 
value taken, and the result displayed. (See Table 4.1) The graph in Figure 
3.^ was then used to classify the MCD's. The appropriate (A), (k), and 
1950 density were then substituted into Equation (2.4) and the predicted 
density calculated for each MCD in 1970, 1980, and 1990. (m = 2,3,4) 

The results of these calculations along with the actual census count for 
the MCD in 1970 are displayed in Table 4.2. 

RESULTS OF THE TEST 

The accuracy of the method was then tested by calculating the per- 
centage error between the 1970 forecast and the actual 1970 oopulation. 

The results of these calculations are displayed in Table 5.1. A summary 
of the magnitude and direction of the errors is displayed in Table 5.2. 

All of the large errors (over 50%) exceot one are under estimates. 

The one over estimate resulted from a change in the boundaries of the 
MCD. The large under estimates are all for MCD's which experienced 
extremely high growth rates between 1950 and 1970. Greenberg et a1 
recognize this weakness and suggest that two to three time oeriods are 
necessary for the model to reach the actual growth of the MCD.*^ A method 
for correcting the error based on proportional replacement of county 
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TABLE 4.1 Parameters used to determine MCD class 


MCD NAME 

LOG (1950 Density) 

LqG ( 50-70 Growth) 

Allen 

3.83 

.070 

Wabash 

3.79 

.021 

Jackson 

4 60 

.063 

Patterson 

3.77 

.045 

Brown 

4.17 

.010 

Wayne 

4.58 

.012 

Washington 

3.54 

-.004 

Greenville Twp. 

5.41 

.016 

Adams 

4.35 

.036 

Liberty 

3.63 

.010 

Neave 

3.75 

.155 

Franklin 

3.92 

-.013 

Harrison 

4.11 

.040 

Butler 

3.67 

.099 

Twin 

4.65 

.109 

Monroe 

4.00 

.093 

Arcanum 

7.31 

.141 

Versailles 

7.41 

.160 

Union City 

8.02 

.056 

Greenville 

7.54 

.181 

Ansonia 

7.75 

.091 

Bath 

6.47 

.378 

Beaver Creek 

4. 68 

1.096 

Ceasar Creek 

3.52 

.082 

Cedarviile Twp. 

4.11 

.199 

Jefferson 

3.69 

.011 

Miami 

5.09 

.178 

New Jasper 

3.62 

.185 

Silver Creek 

4.53 

.110 

Spring Valley 

3.83 

.145 

Sugar Creek 

4.09 

.114 

Xenia Twp. 

5.86 

-.325 

Cedarviile 

7.47 

.341 

Jamestown 

7.85 

.153 

Vellow Springs 

7.37 

.261 

Bellbrook 

4.87 

.689 

Fairborn 

7.01 

.938 

Xenia 

7.64 

.396 

Newberry 

4.93 

.078 

Washington 

6.63 

.089 

Spring Creek 

4.08 

.165 

Newton 

4.15 

.054 

Concord 

6.15 

.157 

Staunton 

3.91 

.659 

Brown 

3.93 

.016 

Lost creek 

3.71 

.088 

Union 

4.85 

.279 

Monroe 

5.15 

.361 

Bethel 

4.22 

.326 

Elizabeth 

3.66 

.120 

Pleasant Hill 

7.60 

.044 


OKlGlN^O. r.v''’ ’ 
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MCD NAME 


West Milton 

Piqua 

Troy 

Tipp City 

Bradford 

Covington 

Butler 

Clay 

German 

Harrison 

Jackson 

Jefferson 

Madison 

Mad River 

Miami 

Perry 

Randolph 

Wayne 

Brookville 

Germantown 

Trotwood 

West Carrolton 

Union 

Centerville 

Englewood 

Dayton 

Kettering 

Miamisburg 

Oakwood 

Jefferson 

Monroe 

Harrison 

Jackson 

Washington 

T..in 

Lanier 

Israel 

Bomers 

Gratis Twp. 

New Paris 
Lewis burg 
Camden 

West Alexandria 
Eaton 


LOG (1950 Density) 


7.18 

8.06 

7.51 
7.11 
7.70 

7.43 
5.33 
4.94 
4.94 
8.05 
4.46 

5.97 

6.20 

7.32 
6.50 
4.56 
5.53 

4.44 

6.98 

6.52 

5.33 

6.21 

5.38 

4.55 
5.27 
8.63 

7.09 
6.83 

8.38 

4.22 

4.09 
4.58 

3.56 
4.97 
4.04 
4.26 
3.62 
4.20 

4.23 
7.35 
7.73 
6.83 
7.81 
6.93 


LOG (50-70 Growth) 


.322 

.090 

.267 

.239 

-.016 

.089 

.709 

.229 

.212 

.213 

.105 

-.029 

.466 

.547 

.446 

.331 

.885 

2.051 

.503 

.281 

1.330 

.662 

1.693 

1.909 

1.843 

-.0006 

.731 

.512 

.021 

.238 

.045 

.114 

.005 

.138 

.151 

.169 

.038 

.130 

.244 

.269 

.123 

.178 

.145 

.190 
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TABLE 4,2 Newling Method Population Projections for Miami Valley Planning Region 


MCD NAME 

1970 

Count 

1970 

Esti/nit'.‘j 

1980 

Estimate 

1990 

Estimate 

DARKE COUNTY 

Allen 

1517 

1580 

1868 

2045 

Wabash 

1022 

1192 

1434 

1584 

Jackson 

2978 

2308 

2035 

1901 

Patterson 

1448 

1622 

1965 

2181 

Brown 

1931 

1970 

2047 

2089 

Wayne 

3736 

3471 

3966 

4234 

Washington 

1074 

1456 

1927 

2244 

Greenville Twp. 

15849 

9606 

7989 

7300 

Adams 

3027 

2728 

2644 

2599 

Liberty 

1278 

1625 

2030 

2378 

Neave 

1311 

1213 

. :83 

1654 

Franklin 

1237 

1461 

1669 

1794 

Harrison 

2205 

2176 

2316 

2396 

Butler 

1623 

1715 

2162 

2452 

Twin 

3623 

3312 

3683 

3885 

Monroe 

1634 

1529 

1699 

1799 

Arcanum 

1993 

1933 

2676 

4204 

Versailles 

2441 

2241 

3010 

4535 

Union City 

1808 

1754 

1956 

2276 

Greenville 

12380 

10656 

13773 

19671 

Ansonia 

1044 

1007 

1220 

1593 

GREENE COUNTY 

Bath 

38474 

30498 

54654 

122903 

Beaver Creek 

26555 

12001 

37079 

177705 

Ceasar Creek 

1071 

1240 

1655 

1936 

Cedarville Twp. 

3346 

2475 

2628 

2715 

Jefferson 

1179 

1459 

1825 

2060 

Miami 

5848 

3936 

3703 

3595 

New JaSper 

1085 

1004 

1287 

1473 

Silver Creek 

2907 

2789 

3249 

3500 

Spring ValUiy 

2136 

2941 

2294 

2511 

Sugar Cree^k 

8276 

4047 

14989 

92387 

Xenia Twp. 

7912 

11977 

8372 

7026 

Cedarville 

2342 

1576 

2076 

3044 

Jamestown 

1790 

1512 

1779 

2229 

Yellow Springs 

4624 

3610 

4904 

7504 

Bellbrook 

1268 

918 

2675 

11816 

Fairborn 

32264 

10539 

16054 

28620 

Xenia 

25373 

15153 

18995 

26002 

MIAMI COUNTY 
Newberry 

6598 

5675 

5673 

5673 

Washington 

22402 

9176 

4776 

3469 

Spring Creek 

2123 

1684 

1810 

1882 

Newton 

2947 

2788 

2922 

2996 

Concord 

19056 

8443 

5275 

4189 

Staunton 

3863 

2081 

3083 

3737 

Brown 

1621 

1799 

2049 

2199 

Lost Creek 

1409 

1493 

1851 

2080 

Union 

9413 

5910 

6081 

6167 

Monroe 

9170 

4472 

4113 

3948 . 

Bethel 

4284 

3260 

4268 

4869 


/" n 



TABLE 4.2 (con't) 



1970 

1970 

1980 

1990 

MCD NAME 

Count 

Estimate 

Estimate 

Estimate 

MIAMI COUNTY (con't) 




Elizabeth 

1456 

1487 

1879 

2134 

Pleasant Hill 

1025 

1114 

1..11 

1960 

West Milton 

3696 

2732 

3934 

6539 

Piqua 

20741 

18690 

20564 

23483 

Troy 

17186 

12913 

16851 

24391 

Tipp City 

5090 

4366 

6430 

11009 

Bradford 

1240 

1487 

1828 

2436 

Covington 

2575 

2673 

3576 

5324 

MONTGOMERY COUNTY 





Butler 

19890 

12673 

32110 

11682' 

Clay 

7438 

4886 

4865 

4856 

German 

7102 

4786 

4760 

4747 

Harrison 

34176 

30046 

33182 

38087 

Jackson 

5823 

3781 

4529 

4947 

Jefferson 

11790 

8040 

5388 

4430 

Madison 

29087 

21440 

41624 

104614 

Mad River 

43881 

22523 

31087 

48638 

Miami 

38705 

27526 

48807 

108148 

Perry 

6020 

3943 

4542 

4868 

Randolph 

20971 

10195 

24331 

81459 

Wayn. 


4561 

15158 

80394 

Brookville 

4403 

2592 

3966 

7160 

Germantown 

4088 

1262 

685 

508 

Trotvraod 

6997 

2034 

5286 

19267 

West Carrolton 

10748 

4628 

8961 

22437 

Union 

3654 

716 

1790 

6393 

Centerville 

10333 

1919 

6180 

31377 

Englewood 

7885 

1342 

3467 

12955 

Dayton 

243601 

230801 

213901 

192244 

Kettering 

69999 

29471 

43683 

75462 

Miami sburg 

14797 

8901 

14294 

27600 

PREBLE COUNTY 





Jefferson 

3598 

2388 

2334 

2461 

Monroe 

2272 

2237 

2394 

2482 

Harrison 

4251 

3018 

2678 

2511 

Jackson 

1211 

1598 

2098 

2431 

Washington 

7748 

5880 

5796 

5755 

Twin 

2612 

2166 

2369 

2488 

Lanier 

3512 

2581 

2592 

2598 

Israel 

1452 

1754 

2253 

2580 

Somers 

2973 

2397 

2465 

2502 

Gratis Twp. 

3782 

2474 

2513 

2534 

New Paris 

1692 

1310 

1792 

2767 

Lewi sburg 

1553 

1413 

1728 

2275 

Camden 

1507 

1525 

2449 

4729 

VJest Alexandria 

1553 

1340 

1594 

2027 

Eaton 

6020 

5834 

9081 

16791 


0 . ' 
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TABLE 5.1 Results of the Newling Model for Miami Valley Planning Region 


MCD NAME 

Class 1950 

1970 

1970 

% Error 


Density 

Pop. 

Est. 



Allen 

1 

46.3 

1517 

1580 

4 

Wabash 

1 

44.2 

1022 

1192 

17 

Jackson 

1 

99.1 

2978 

2308 

-23 

Patterson 

1 

43.4 

1448 

1622 

12 

Brovm 

1 

64.7 

1931 

1970 

2 

Wayne 

2 

97.7 

3736 

3471 

- 7 

Washington 

1 

34.5 

1075 

1456 

36 

Greenville Twp. 

2 

223.1 

15849 

9606 

-39 

Adams 

1 

77.5 

3017 

2728 

-10 

Liberty 

1 

37.5 

1278 

1625 

27 

Weave 

1 

42.4 

1311 

1213 

- 7 

Franklin 

1 

50.6 

1237 

1461 

18 

Harrison 

1 

42.4 

2205 

2167 

- 1 

Butler 

1 

39.1 

1623 

1715 

6 

Twin 

2 

104.7 

3623 

3312 

- 9 

Monroe 

1 

54.3 

1634 

1529 

- 6 

Arcanum 

3 

1499.5 

1993 

1933 

- 3 

Versailles 

3 

1654.3 

2441 

2241 

- 8 

Union City 

3 ■ 

3044 . 2 

1808 

1754 

- 3 

Greenville 

3 

1876.8 

12380 

' 10655 

-14 

Ansonia 

3 

2319.3 

1044 

1007 

- 4 

Bath 

3 

643.4 

38474 

30498 

-21 

Beaver Creek 

3 

108.0 

26555 

12001 

-55 

Ceasar Creek 

1 

33.7 

1071 

1240 

16 

Cedarville Twp. 

1 

61.6 

3346 

2475 

-26 

Jefferson 

1 

40.0 

1179 

1459 ' 

24 

Miami 

2 

161.9 

5848 

3936 

-33 

Wew Jasper 

1 

37.4 

1085 

1004 

- 7 

Silver Creek 

2 

93.0 

2907 

2789 

- 4 

Spring Valley 

1 

46.3 

2136 

1941 

- 4 

Sugar Creek 

3 

59.7 

8276 

4047 

-51 

Xenia Twp. 

2 

350.4 

7912 

11977 

51 

Cedarville 

3 

1763.1 

2352 

1576 

-33 

Jamestown 

3 

2554.3 

1790 

1512 

-16 

Yellow Springs 

3 

1595.0 

4624 

3610 

-22 

Bellbrook 

3 

130,9 

1268 

918 

-28 

Fairborn 

3 

1112.6 

32264 

10589 

-67 

Xenia 

3 

2074.3 

25373 

15153 

-40 

Newberry 

2 

138.4 

6598 

5675 

-14 

Washington 

2 

753.9 

22402 

9176 

-59 

Spring Creek 

1 

59.2 

2123 

1684 

-21 

Newton 

T 

i. 

63.3 

2947 

2788 

- 5 

Concord 

2 

469.7 

19056 

8448 

-56 

Staunton 

2 

49.8 

3863 

2081 

-46 

Brown 

1 

51.0 

1621 

1700 

11 

Lost Creek 

1 

40.8 

1409 

1493 

6 

Union 

2 

128.3 

9413 

5910 

-37 
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! TABLE 5.1 (con't) Results of the Mewling Itodel for Miami Valley Planning Region 


MCD NAME 

Class 1950 

1970 

1970 

% Error 


Density 

Pop. 

Est. 



Monroe 

2 

171.8 

9170 

4472 

-51 

Bethel 

2 

68.7 

4284 

3260 

r24 

Elizabeth 

1 

38.8 

1456 

1487 

2 

Pleasant Hill 

3 

2005. 3 

1025 

1114 

9 

West Milton 

3 

1317.0 

3696 

2732 

“26 

Piqua 

3 

3180.3 

20741 

18690 

-10 

Troy 

3 

1818.5 

17186 

12913 

-25 

Tipp City 

3 

1223.7 

5090 

4366 

-14 

Bradford 

3 

2209.8 

1240 

1487 

20 

Covington 

3 

1691.1 

2575 

2673 

4 

Butler 

3 

206.9 

19890 

12673 

-36 

Clay 

2 

139.7 

7438 

4886 

-34 

German 

2 

140.2 

7102 

4786 

-33 

Harrison 

3 

3142.0 

34176 

30046 

-12 

Jackson 

2 

86.5 

5823 

3781 

-35 

Jefferson 

2 

390.9 

11790 

8040 

-32 

Madison 

3 

495.0 

29087 

21440 

-26 

Mad River 

3 

1513.6 

43881 

22253 

-49 

Miami 

3 

666.2 

38705 

27256 

-29 

Perry 

2 

95.7 

6020 

3943 

-34 

Randolph 

3 

251.7 

20971 

10195 

-51 

Wayne 

3 

85.1 

27975 

4561 

-84 

Brookville 

3 

1079.7 

20971 

10195 

-51 

Germantown 

2 

676.4 ' 

4088 

1262 

-69 

Trotwood 

3 

206.1 

6997 

2084 

-70 

West Carrolton 

3 

499.4 

10478 

4628 

-57 

Union 

3 

216.1 

3644 

716 

-80 

Centerville 

3 

94.4 

10333 

1919 

-81 

Englewood 

3 

194.3 

7885 

1342 

-83 

Dayton 

3 

5588.2 

243601 

230801 

- 5 

Kettering 

3 

1198.3 

69999 

29571 

-58 

Miaraisbunj 

3 

921.6 

14797 

8901 

-40 

Oakwood 

3 

4364.7 

10095 

9683 

- 4 

Jefferson 

1 

67.9 

3598 

2388 

-34 

Monroe 

1 

59.9 

2272 

2237 

- 2 

Harrison. 

1 

97.5 

4251 

3018 

-29 

Jackson 

1 

35.2 

1211 

1598 

32 

Washington 

2 

143.5 

7748 

5880 

-24 

Twin 

1 

56.6 

2612 

2166 

-17 

Lanier 

1 

70.7 

3512 

2581 

-27 

Israel 

1 

37.3 

1452 

1754 

21 

Somers 

1 

66.5 

2973 

2397 

-19 

Gratis Twj'. 

1 

63.6 

3782 

2474 

-35 

New Paris 

3 

1560.5 

1692 

1310 

-23 

Lewisburg 

3 

2275.1 

1553 

1418 

- 9 

Camden 

3 

921.3 

1507 

1525 

1 

West Alexandria 

3 

2466.2 

1553 

1340 

-14 

Eaton 

3 

1020.6 

6020 

5834 

- 3 
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total errors In presented is Equation (5.1). 

■ 1970 MCD Count 

Eq. (5.1) CE = 1970 County Count *(Total County Error); 

If FE>1970 MCD Count. CE = FE - CE' ; 

If FE<1970 MCD Count, CE = FE + CE' . 

where CE* is the corrected estimate. 

1970 MCD Count is the census count for the MCD in 1970, 

1970 County Count is the census count for the county in 1970, 

Total Error is the difference between the sum of the MCD estimates 

for the county and the 1970 County Count, and 

FE is the original estimate of the MCD population. 

This method can be extended for any decade for which a reliable 
county total estimate exists by replacing 1970 MCD Count with the Newling 
model estimate and 1970 County Count with the county total estimate. 

The value of the Newling model lies in its simplicity of execution 
once the parameters have , been determined. It is not a rigorous forecasting 
method, but it serves the "quick and dirty" puroose well. 



FOOTNOTES 


^Michael R. Greenberg, Kreuckenberg, Donald A., and Mautner, 

Richard, Long-Range Population Pro.iectlons for Minor Civil Divisions . 
Center for Urban Policy Research, Rutgers' University, New Brunswick, 
N. J., 1973. 

^Greenberg et al, pp. 16-18. 

^End of direct quotation 
^Greenberg et al, P. 32. 
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